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ABSTRACT 

The first-year WMAP data, in combination with any one of a number of other 
cosmic probes, show that we live in a flat A-dominated CDM universe with m ~ 0.27 
and J7a ~ 0.73. In this model the late-time action of the dark energy, through the 
integrated Sachs- Wolfe effect, should produce CMB anisotropies correlated with matter 
density fluctuations at z < 2 (Crittenden & Turok 1996). The measurement of such a 
signal is an important independent check of the model. We cross-correlate the NRAO 
VLA Sky Survey radio source catalog (Condon et al. 1998) with the WMAP data in 
search of this signal, and see indications of the expected correlation. Assuming a flat 
ACDM cosmology, we find J7a > (95% CL, statistical errors only) with the peak of 
the likelihood at fl\ = 0.68, consistent with the preferred WMAP value. A closed 
model with Q m = 1.28, h = 0.33, and no dark energy component (Q\ = 0), marginally 
consistent with the WMAP CMB TT angular power spectrum, would produce an anti- 
correlation between the matter distribution and the CMB. Our analysis of the cross- 
correlation of the WMAP data with the NVSS catalog rejects this cosmology at the 3a 
level. 
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1. 



Introduction 



The recent WMAP results (Bennett et al. 2003c) place on a firm foundation the emerging 
standard model of cosmology: a flat adiabatic A-dominated CDM universe. However, the deficit 
of power in the CMB anisotropy spectrum on large angular scales (Bennett et al. 2003b; Hinshaw 
et al. 2003; Spergel et al. 2003) is surprising given that ACDM models predict an enhancement 
at I < 10 due to the late-time integrated Sachs- Wolfe (ISW) effect (Sachs & Wolfe 1967). The 
likelihood of observing so little power due to sample variance is only 0.15% (Spergel et al. 2003). 
Thus new physics may be indicated since the nature of the dark energy is poorly understood. Cross- 
correlating the CMB with radio sources provides a direct test for the recent acceleration predicted 
by the favored ACDM model and observed by the Type la supernovae experiments (Perlmutter 
et al. 1999; Riess et al. 1998). Thus, this test is an important check of the standard model. 

Any recent acceleration of the universe will cause local gravitational potentials to decay. This 
decay is then imprinted on the CMB as the photons, which were blue-shifted on infall, suffer less 
of a red-shift as they climb out of the potential well. This produces temperature perturbations 



where $ is the Newtonian gravitational potential, rj is conformal lookback time, and the integral 
runs from today (rj = 0) to the CMB decoupling surface at z^ cc = 1089. Figure 1 shows the recent 
evolution of <1> for a variety of cosmological models. Since $ is related to the matter distribution via 
the Poisson equation, tracers of the mass will be correlated with the CMB through the late-ISW 
effect. In this paper we correlate the NRAO VLA Sky Survey (NVSS) source catalog (Condon et al. 
1998) with the WMAP CMB map in search of the late-ISW effect. Boughn & Crittenden (2002) 
performed a similar analysis using the COBE/DMR map, and found no correlation. However, 
a recent re-analysis by the same authors using the WMAP data did see evidence for a correla- 
tion between the CMB, NVSS, and the hard X-ray background observed by the HEAO-1 satellite 
(Boughn & Crittenden 2003). Here we focus on the implications of this correlation for dark energy, 
specifically Oa- 

We take as our fiducial model the best fit power-law ACDM model to the combined WMAP, 
CBI, ACBAR, 2dFGRS, and Lya data sets, with values uj m = 0.133, oj b = 0.0226, n s = 0.96, 
h = 0.72, A = 0.75, and r = 0.117 (Spergel et al. 2003, Table 7). Cosmological parameters are 
drawn from this set unless otherwise noted. 

In §2 we briefly describe the NVSS source catalog and its auto-correlation function. In §3 we 
describe the cross-correlation between the WMAP CMB map and the NVSS source map, and relate 
it to the late-ISW effect. We conclude in §4 and discuss effects which could mimic the observed 
signal. 
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2. The NVSS Source Catalog 

The NRAO VLA Sky Survey (NVSS) is a 1.4 GHz continuum survey, covering the 82% of 
the sky with 5 > —40° (Condon et al. 1998). The source catalog contains over 1.8 x 10 6 sources, 
and is 50% complete at 2.5 mJy Nearly all the sources away from the Galactic plane (\b\ > 2°) 
are extragalactic. The bright sources are predominantly AGN-powered radio galaxies and quasars, 
while at weaker fluxes the sources are increasingly nearby star-forming galaxies. 

Galaxy counts are a potentially biased tracer of the underlying matter distribution, and thus 
the projected number density of NVSS sources per steradian, n(z,h), is related to the matter 
distribution 5(z,h) = 5p/p via 

n(z, n) dz dCl = (1 + b r (z)8(z, n)) dz dCl (2) 

where dN/dz is the mean number of sources per steradian at a redshift z and b r (z) is the radio 
galaxy bias parameter. Thus the observed fluctuation on the sky in projected source counts is given 
by 

SN(h) = j ' dzb r {z)^5{z,h). (3) 

Since we are only interested in clustering on large scales, and hence the linear regime, the evolution 
of 5 factors as 5(k, z) = D(z)5(k), where S(k, z) is the Fourier transform of the matter distribution, 
5(h) = S(k,0) its current value, and D(z) is the linear growth factor (Peebles 1993, 5.111). While 
generally a function of time, we take the bias to be constant, as the determination of dN/dz is 
uncertain and we only consider a modest redshift range (0 < z < 2). 

While the individual NVSS source redshifts are unknown, for our purposes we need only the 
overall redshift distribution dN/dz for the NVSS. We adopt the favored model of Dunlop & Peacock 
(1990) (model 1, MEAN-z data), which was found by Boughn & Crittenden (2002) to best reproduce 
the NVSS auto-correlation function. Like most models, it divides the sources by spectral index a 
(flux S oc v~ a ) into two populations, flat-spectrum (a ~ 0) and steep-spectrum (a ~ 0.8) sources. 
The model is plotted in Figure 2. We limit the model to 0.01 < z < 5; the lower limit corresponds 
to a distance of « 42Mpc. The small peak at z ~ 0.05 is spurious (due to a breakdown in the 
DP90 fitting function), but has only a minor effect on the integrated predictions. 

Given a cosmology, we can determine the radio bias b r from the amplitude of the NVSS auto- 
correlation function (ACF), by comparing the ACF to the unbiased prediction 

C NN (9) = (5N(n)5N(n>)) = £ ^[ftff c/™ m ■ n>) (4) 

where [bf] 2 is the pixel window function. Substituting equation (3) into this expression and Fourier 
transforming, 

cr=^Jf^(k)[f?(k)] 2 (5) 
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where A$(k) = k 3 Ps(k)/2ir 2 is the logarithmic matter power spectrum and Ps(k) = (\S(k)\ 2 ). We 
use the WMAP normalization A^(fc ) = 2.95 x W~ 9 A where k = 0.05 Mpc -1 (Verde et al. 2003), 
giving us 5h = 6.1 x 10 -5 for our fiducial model. The filter function is given by 

f?{k)=b r j dz^Diz^kv). (6) 

where ji{x) is the spherical Bessel function. Note that the bias we measure is complicated by the 
uncertainty in dN/dz; errors in the normalization are absorbed into the bias. 

To calculate the observed NVSS auto-correlation function (ACF), we made a HEALPix (Gorski 
et al. 1998) resolution-5 map of 5N(h), which has 12 288 1.8° square pixels 10 . As a precaution, we 
removed the 3 x 10 5 sources from the catalog which were resolved. The mean source count per 
pixel is 147.9, leading to a Poisson uncertainty of « 8%. The ACF is estimated as 

C™y k ) = £ N iNj v?uf/ £ wf W f (7) 

where Ni is the number of sources in pixel i and the sum is over all pixel pairs separated by 
9 k - A6/2 < 9 < 9 k + A6/2. The bin width A6> is 2°. We mask out pixels at low Galactic latitude 
< 10°) and those unobserved by the survey (5 < —37°); the weights wf are determined by 
the mask. The NVSS ACF is shown is Figure 2. The 8 = 0° bin is corrected for Poisson noise by 
subtracting the mean number of sources per pixel. For the fiducial model parameters (which are 
used to calculate D[z) and r](z)), the derived bias is 1.7. This is somewhat higher than the value 
of 1.6 found by Boughn <fc Crittenden (2002). However, they assumed a scale- invariant spectrum 
(i.e., n s = 1), and changing n s by ±0.03 changes the bias by =f0.05. 

As noticed by Boughn & Crittenden (2002), the NVSS catalog mean source density varies with 
declination, introducing a spurious signal into the auto-correlation function. They corrected for 
this by adding and subtracting random sources from the map until the structure was removed. We 
considered two simple corrections; both broke the sources into sin(<5) strips of width 0.1. The first 
method subtracted the mean from each strip; the second scaled each strip by the ratio of the global 
mean to the strip mean. Since the corrections are small, both produced similar results (Fig. 3). 



3. WMAP— NVSS Cross-Correlation 

Combining equation (1) and equation (3) we can calculate the expected cross-correlation spec- 
trum between the NVSS catalog and the CMB: 

Cr = <«> = 4vr I f Al(k) fl N (k)f?(k) (8) 



10 For more information on HEALPix visit http://www.eso.org/science/healpix/. The resolution of a HEALPix 
map indicates its pixel count. A resolution-r map has 12iV s 2 ide pixels, where iV s idc = 2 r . A resolution-(r + 1) map is 
created by dividing each resolution-r pixel into four subpixels. 
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where f t N and ff are the NVSS and ISW filter functions. The ISW filter function is derived 
analogously to the NVSS filter function. The local gravitational potential is related to the matter 
distribution vici th.6 Poisson equation — AttGqP 1 Pm^mi 

where the gradient is taken with respect 

to comoving coordinates. Fourier transforming, we have 

= -\n m {H Q lk?g{rj)5{k) (9) 

where Hq is the Hubble constant, £l m HQ = 87rGp^/3, and 17(77) = D(rj)/a(rj) is the linear growth 
suppression factor. Thus 

ff{k) = m m {H /kf j d V ^ji(kv)- (10) 

We use the fitting function for g{rj) provided by Carroll et al. (1992). In a flat Q m = 1 universe, 
g{rf) is constant, and thus there is no ISW effect and hence no correlation between the CMB and the 
local matter distribution. In ACDM universes, D{rj) approaches a constant during A-domination, 
leading to a decay of 17(77) as time increases. 

We computed the WMAP-N VSS cross-correlation function (CCF) as 

C NT (0 k ) = £ N^wfwj/ £ wfwj (11) 

where Tj is the CMB map, Nj the NVSS map (described in the previous section), and the sums are 
over all pixel pairs separated by k - AO/2 < < k + AO/2. The bin width AO is 2°. As before, 
we work at HEALPix resolution-5. Since we are working at large scales where the detector noise 
is negligible, we use the WMAP internal linear combination (ILC) map for the CMB (Bennett 
et al. 2003b). We limit residual foreground contamination by masking the map with the WMAP 
KpO galaxy mask and the WMAP source mask (Bennett et al. 2003a). The CMB weight wj is 
the number of unmasked resolution-9 subpixels of the resolution-5 pixel i. The CCF is plotted in 
Figure 3. The CCF is insensitive to the form of the NVSS declination correction. 

Assessing the significance of the CCF is complicated by the high degree of correlation between 
points. Accidental alignments between the NVSS map and the CMB fluctuations at the decoupling 
surface (which are uncorrelated with those generated by the late-ISW) can produce spurious cor- 
relations. We quantified this uncertainty with Monte Carlo simulations, creating 500 realizations 
of the CMB sky drawn from the Cf T power spectrum for our fiducial model with our estimate of 
the NVSS-correlated late-ISW contribution subtracted. The covariance matrix S was calculated 
from the resulting CCFs, keeping the NVSS map fixed. Table 1 shows the 0° < < 20° subma- 
trix of S, showing that the CCF points are highly correlated. We define x 2 = 5C T T I ~ 1 5C where 
5C = C NT — C NT is the difference between observed (eq. [11]) and model (eq. [8]) correlation 
functions, and we limit < 20°. For the null model of no correlation (C NT = 0), Xo = 17.2. Since 
there are 10 degrees of freedom this a 1.8<r deviation. 

How does the measured CCF constrain Rather than exploring the full parameter space, 
we assume a flat universe with fixed ujb and explore the locus of values of Q. m and h consistent with 
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the measured location of the first acoustic peak of the CMB TT anisotropy power spectrum. The 
first peak position is set by the angular scale of the sound horizon at decoupling, 9 a, which Page 
et al. (2003) found to be 9 A = 0.6°. Percival et al. (2002) showed that 9 A » 0.85°^ 14 /i a48 ; this 
is the horizon angle degeneracy. The normalization A is varied to fix the amplitude of the first 
peak; from a fit to CMBFAST (Seljak & Zaldarriaga 1996) spectra we found A oc Q^ 248 along the 
horizon degeneracy. The results are shown in Figure 4. The difference in \ 2 between models with 
Qa = (i.e., no correlation) and = 0.68 (the minimum, with Xmin = 12.5) is A% 2 = 4.7. Since 
we are varying a single parameter, the significance is y A% 2 ; thus Q\ > is preferred at the 2.2<r 
level. 

The WMAP team imposed a prior on the Hubble constant, h > 0.5, in determining the 
cosmological parameters (Spergel et al. 2003). While lower values of the Hubble constant would 
contradict a host of other experiments, especially the Hubble Key Project (Freedman et al. 2001), 
models with very low h and fi m ss 1.3 are marginally consistent (A% 2 = 4.9) with the WMAP TT 
and TE angular power spectra. Since these universes are closed and matter dominated the growth 
factor D(a) grows faster than a, and g(a) is better termed the linear growth enhancement factor. 
Thus we would expect to observe an anti-correlation between the NVSS and CMB maps, since 
dg/dr] in equation (10) changes sign (see Figure 1). For O m = 1.28 and h = 0.33, % 2 = 24.2 (with 
a bias of 2.7). Based on the cross-correlation analysis alone, this model is disfavored by more than 
3d. 

4. Discussion 

The recent acceleration of the universe due to dark energy should correlate large-scale CMB 
anisotropies with fluctuations in the local matter density through the late-time integrated Sachs- 
Wolfe effect. We have correlated the NVSS radio source catalog with the CMB anisotropies observed 
by the WMAP satellite, and find that f^A > is preferred at the 95% confidence level (A% 2 = 4.7), 
considering statistical errors only. The statistical uncertainty is due to accidental alignments with 
the background primary anisotropies generated at decoupling. The likelihood peaks at Qa = 0.68, 
consistent with the value derived from the CMB angular power spectrum. 

The correlation between the NVSS source count and WMAP CMB maps appears robust. We 
interpret it as arising from the late-ISW, but other effects could correlate the two maps. For 
instance, obscuration by dust clouds tends to reduce the number of sources observed in their 
direction. We cross-correlated the NVSS map with the E(B — V) extinction map of Schlegel et al. 
(1998) and see evidence for a small negative correlation at separations 9 < 20°. However, since 
the extinction map is positively correlated with the CMB map due to dust emission, this effect 
has the wrong sign to mimic the late-ISW. The extinction correction is estimated as te{9) = 
(EN)(ET) /(EE), where N, T, E are the NVSS, CMB, and extinction maps and (XY) denotes the 
correlation between maps X and Y evaluated at separation 9. We find a value of rg(0) « —4 jjK cnts 
(compared to 22 /iKcnts for the CCF) and the correction is negligible for 9 > 15°, except for a few 
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glitches when (EE) crosses zero. Subtracting te from the CCF, we find the preferred value of f^A 
increases to 0.76; at the minimum Xmin = 14.8, and Ax 2 = 6.1. 

A potentially more serious concern is that the correlation is due to microwave emission by the 
sources themselves. However, if this were the case then the CCF should have a similar angular 
profile as the auto-correlation function (ACF). Yet while the ACF falls steeply with increasing 
separation [C NN (0°)/C NN (3°) ~ 5], the CCF does not [C NT {0°)/C NT (3°) ~ 1]. The lack of an 
enhanced signal in the zero lag CCF bin thus argues against any significant microwave emission 
from the NVSS radio sources. 

We are indebted to Steve Boughn for significant and helpful discussions throughout the prepa- 
ration of this paper. We also thank Hiranya Peiris and Licia Verde for useful comments. The 
WMAP mission is made possible by the support of the Office of Space Sciences at NASA Head- 
quarters and by the hard and capable work of scores of scientists, engineers, technicians, machinists, 
data analysts, budget analysts, managers, administrative staff, and reviewers. 
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Table 1. WMAP-NVSS CCF Correlation Matrix 



Bin 


1° 


3° 


5° 


7° 


9° 


11° 


13° 


15° 


17° 


19° 


1° 


183 


159 


141 


131 


115 


98.7 


87 


73.7 


58.7 


49.7 


3° 




161 


140 


136 


119 


104 


94.3 


82.8 


67.6 


59.6 


5° 






135 


128 


113 


97.8 


86.4 


74.8 


60.4 


52.5 


7° 








135 


120 


106 


96.7 


86.4 


71.9 


64.2 


9° 










113 


102 


92.6 


83.7 


70.3 


63.4 


11° 












95.9 


89.2 


81.6 


69.6 


63.5 


13° 














88.3 


82.7 


71.7 


66.3 


15° 
















81.7 


72.3 


67.6 


17° 


















67.5 


64.2 


19° 




















64.6 



Note. - ■ The correlation matrix S of the WMAP-NVSS cross-correlation function due to accidental 
alignments with the anisotropies produced at decoupling. Units are (/zKcnts) 2 . 
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CI =0.27, A = 0.73 




Q = 1 .28, A = 

m 



0.0 0.5 1.0 1.5 2.0 

Redshift 



Fig. 1. — The gravitational potential f as a function of redshift z for a variety of cosmological 
models. The models are normalized to unity at z = 0. 
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Fig. 2. — The adopted dN/dz model (RLF1) for the distribution of NVSS sources from Dunlop & 
Peacock (1990, DP90), normalized to integrate to unity (left panel). The small blip at z ~ 0.05 is 
spurious, due to breakdown in the DP90 fitting function. Also plotted is the "luminosity /density- 
evolution" (LDE) model also from DP90, which is a poor fit to the observed auto-correlation 
function (right panel). 
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Fig. 3.— The WMAP-NVSS cross-correlation function (CCF). The CCF is insensitive to the 
details of the declination correction. Two simple methods are compared; both broke the sources 
into sin((5) strips of width 0.1. The first (diamonds) subtracted the mean from each strip. The 
second (triangles) scaled each strip by the ratio of the global mean to the strip mean. The cross 
points are uncorrected, showing the correction is only important for 6 > 25. We used the WMAP 
internal linear combination (ILC) CMB map; substituting the map of Tegmark et al. (2003) instead 
produces the same results (square points). The solid and dashed lines are derived from the diagonal 
elements of the correlation matrix due to accidental alignments; they would be the la and 2a 
contours in the absence of off-diagonal correlations. The points, however, are highly correlated as 
shown in Table 1. 
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Fig. 4. — Effect of varying S7a on the cross-correlation function. In all panels we assume a flat 
universe with fixed uib, and trade off between Oa and h by keeping the combination Q, m h? A constant; 
when = 0, h = 0.48. Panel (a) shows the inferred radio bias as a function of Q\. Panel (b) shows 
the bias-corrected NVSS auto-correlation function (ACF) compared with the measured ACF. Panel 
(c) shows the predicted cross-correlation function (CCF) for a range of values of Qa, compared with 
the measured CCF (diamonds). The amplitude of the predicted CCF is proportional to Qa, which 
is stepped in increments of 0.1 from 0.0 to 0.9. Panel (d) shows the x 2 °f the model CCF as a 
function of Qa- The x 2 was computed using the first 10 points of the CCF (0° < 8 < 20°) and 
Table 1. The minimum Xmin ^ s 

12.5 at n A = 0.68; at S7 A = 0, Xo = 17.2. The dashed line is the 
likelihood oc exp(-x 2 / 2 )- The lcr limits are °- 42 < < 0.86. 



